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Abstract —The theoretical analysis of downlink massive MIMO 
usually assumes the ideal Gaussian channel matrix with asymp¬ 
totic orthogonality of channel vectors associated with different 
users, since it can provide the favorable propagation condition. 
Meanwhile, recent experiments have shown that massive MIMO 
channels between a certain user and massive base station an¬ 
tennas appear the spatial common sparsity (SCS) in both the 
time domain and angle domain. This motivates us to investigate 
whether realistic sparse massive MIMO channels could provide 
the favorable propagation condition, and reveal the capacity gap 
between massive MIMO systems over realistic sparse channels 
and that under the ideal Gaussian channel matrix assumption. 
This paper theoretically proves that channel vectors associated 
with different users in massive MIMO over realistic sparse 
channels satisfy the asymptotic orthogonality, which indicates 
that the favorable propagation condition can also be provided. 
Moreover, the simulation results confirm the theoretical analysis. 

Index Terms —Massive MIMO, favorable propagation condi¬ 
tion, spatial common sparsity (SCS), asymptotic orthogonality. 


I. Introduction 

The theoretical performance of massive MIMO relies on the 
favorable propagation condition, which is equivalent to the 
asymptotic orthogonality of channel vectors associated with 
different users m. For theoretical analysis, the channel is 
usually modeled as the ideal Gaussian channel matrix with 
elements following the mutually independent Gaussian distri¬ 
bution, so that channel vectors associated with different users 
are asymptotically orthogonal and can provide the favorable 
propagation condition for massive MIMO Cl, lEl. 

Recent researches have shown that, due to the finite number 
of significant scatterers and the compact antenna array at the 
base station (BS), realistic massive MIMO channels between a 
certain user and massive base station (BS) antennas exhibit the 
intrinsically spatial common sparsity (SCS) E-ii. However, 
in a-in, the proposed channel estimation schemes only 
consider the massive MIMO channels with SCS in the time 
domain, while Q and El only consider the SCS of massive 
MIMO channels in the angle domain. For massive MIMO 
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channels, both the path delay in the time domain and angle of 
departure (AoD)/angle of arrival (AoA) in the angle domain 
are essential parameters to represent the channels. Therefore, 
in this paper, we will jointly consider the SCS of massive 
MIMO channels in both the time domain and angle domain. 
Clearly, the actual massive MIMO channel experienced by 
practical systems is different from the ideal Gaussian channel 
matrix for theoretical analysis. Such difference motivates us 
to investigate the propagation condition of realistic massive 
MIMO channels with SCS, and reveal the capacity gap be¬ 
tween massive MIMO systems over actual sparse channels 
and that under the ideal Gaussian channel matrix assumption. 
Specifically, we first introduce the SCS of massive MIMO 
channels in both the time domain and angle domain, and 
illustrate the connection between the sparse massive MIMO 
channel model and the frequency-domain channel matrix 
model in massive MIMO-OFDM systems. Then, we prove 
that channel vectors over sparse channels are asymptotically 
orthogonal, where our analysis is based on the widely used 
isotropic uniform linear array (ULA) without loss of generality 
a, Q, 0. The asymptotic orthogonality of different channel 
vectors indicates that massive MIMO channels with SCS could 
provide the favorable propagation condition, which is also 
confirmed by simulations. 

II. System Model 

In typical massive MIMO systems, the BS employs M 
antennas to simultaneously serve K single-antenna users. 
Usually, M is large, and much larger than K, e.g., M = 128 
and iT = 16 H]. In the downlink, the BS transmits data 
X £ to K single-antenna users, and the received signal 

y £ for K users can be expressed as 

y = v7dGx-f n = + n, (1) 

where pd is the transmit power, G £ is the downlink 

channel matrix, n £ is the additive white Gaussian 

noise (AWGN) vector, the diagonal matrix 
with the fcth diagonal element represents the large-scale 
channel fading due to shadow fading and path loss, while H £ 
(^KxM (jenotes the small-scale channel fading matrix due to 
channel variation. Note that when the widely used OFDM is 
considered, O is valid for every frequency-domain subcarrier. 

To fully exploit advantages of massive MIMO, the BS 
requires the downlink channel information for precoding, 
user scheduling, etc. The propagation condition of downlink 
channel matrix determines the achievable capacity of downlink 
massive MIMO, which can be expressed as m 

G = log2det {iK+PdGG*), 


( 2 ) 
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where Ix € is the unit matrix, and (•)* is the 

Hermitian transpose operation. It can be observed from (|2]i 
that the achievable capacity heavily depends on the property 
of channel matrix G. 

For the ideal Gaussian channel matrix, G provides the 
favorable propagation condition, where elements of the small- 
scale channel fading matrix H follow the independent and 
identically distributed (i.i.d.) circularly symmetric complex 
Gaussian distribution CJ\f{Q, 1). It has been shown that row 
vectors of G are asymptotically orthogonal fT), i.e., 

lim GG*/M= lim /M = T>. (3) 

M—¥oo M—¥oo 

Under such ideal Gaussian channel matrix assumption, the 
available capacity has the simple asymptotic form as follows 

K 

log 2 det {Ik + PdGG*) ^ log 2 (1 + pdMp,). (4) 

Such favorable propagation condition can ensure that the inter¬ 
user interference vanishes when M ^ oo m. Accordingly, 
the superiority of massive MIMO in spectral efficiency can be 
guaranteed. For practical massive MIMO with a finite M, the 
favorable propagation condition indicates that the condition 
number of G should be as small as possible IT]. 

III. Asymptotic Orthogonality Analysis for 
Massive MIMO Channels with SCS 

We first illustrate the realistic massive MIMO channels 
with SCS in both the time domain and angle domain. Then 
we theoretically prove that channel vectors are asymptotically 
orthogonal, where the first and second-order moments of the 
inner product of different channel vectors are derived. 

A. Massive MIMO Channels with SCS 

Due to the finite number of significant scatterers, for the 
time-domain channel impulse response (CIR), the majority of 
channel gains are concentrated on a small number of channel 
taps although the delay spread can be large in typical wireless 
communication systems 0. The sparse downlink time-domain 
CIR between the fcth user and the mth BS antenna can be 
expressed as 0 

hkm — ] ,^km ^ Skm Tmaxfs, (5) 

^ - ^Skm = l 

where and are the Sfemth path gain and path delay, 
respectively, Skm is the number of resolvable multipaths, Tmax 
is the maximum delay spread, fs is the system bandwidth, and 
Anax/s is the normalized maximum delay spread. 

Moreover, there are measurements showing that time- 
domain CIRs between different co-located BS antennas and a 
certain user share very similar path delays H), i.e., Skm = Sk 
and Skm = Sk, for 0 < m < M — 1. This is because the scale 
of the compact antenna array at the BS is small compared with 
the large distance of signal transmission, which leads to the 
fact that channels of different transmit-receive antenna pairs 
share very similar scatterers. The inherent sparsity of time- 
domain massive MIMO channels discussed above is referred 
to as the SCS in the time domain. 


Furthermore, for a certain multipath from one scatterer, the 
AoD seen from M antennas at the BS are very similar, which 
is referred to as the SCS in the angle domain of massive 
MIMO channels 0, 0, 0. Due to such SCS in the angle 
domain, the path gains associated with M antennas via the 
Sfeth scatterer can be expressed as the steering vector, i.e.. 


1 


Sjc 


1, 


■Sfc(M-l) 
Kd sin(0f ^ ) 


27r2£isin(<9^«) 


, 27 r(M-l)(isin( 0 ,«) 


,e-^ 


(6) 


where the isotropic ULA with the antenna spacing d is 
considered without loss of generality, 0 :^^“ = is the path 
gain of the first (m = 0) BS antenna, 0^'“ is the AoD of the 
Sfcth propagation path, and A denotes the wavelength. 

The channel model with SCS in both the time domain and 
angle domain is different from the ideal Gaussian channel ma¬ 
trix with the inherent favorable propagation condition widely 
used for theoretical analysis. This motivates us to investigate 
whether the realistic sparse massive MIMO channels meet the 
favorable propagation condition, and reveal the capacity gap 
between massive MIMO over realistic sparse channels and that 
under the ideal Gaussian channel matrix assumption. 


B. Asymptotic Orthogonality Analysis 

According to 0, the channel frequency response over the 
nth subcarrier, denoted by gkmln], can be expressed as 

9km[n] = “fcm exp 

-Sfcm —1 /7\ 

Skm 

= E al-m exp i-j2TrjnTk'^), 1 < n < A^, 

where N is the size of the OFDM symbol and 7 „ = nfs/N. 

Eq. 0 builds the connection between time-domain massive 
MIMO channels and frequency-domain channel matrix. For 
multi-user MIMO-OFDM systems, gkm[n] is the element of 
the fcth row and the mth column in the channel matrix 
G[n] over the nth subcarrier (see G in Eq. 0). Here 
G[n] = [(gi[n])'^,(g 2 [n])'^,--- , (g/f[n])"^]^ and gfc[n] = 
[9ki[n], gk2[n], ■ ■ ■ ,gkM[n]]. Based on 0, we will investigate 
the asymptotic orthogonality of row vectors of the frequency- 
domain channel matrix over sparse massive MIMO channels. 

Due to the SCS in the angle domain as shown in 0, for 
the Sfcth propagation path, the path gain associated with the 
mth BS antenna can be expressed as 

0 < m < M - 1. (8) 

Eurthermore, due to the SCS in the time domain, under the 
far field assumption, the corresponding path delay of the mth 
BS antenna via the Sfcth scatterer can be expressed as 0 


-f mdsin(0^'“)/c, 0 < m < M - 1, (9) 

where is the path delay of the first (m = 0) BS antenna, 
and c is the velocity of electromagnetic waves. 

Consequently, by substituting 0 and 0 into 0, we have 

Sk mdsinfsf^) gj 

gkmln] = E 

, nn') 

Sk mdsin( 6 ,^) ^ ' 

«fc=l 
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where . 

For massive MIMO systems, since the distance between 
different users simultaneously communicating with the BS 
is usually large, their respective CIRs are usually mutually 
independent fl]. Meanwhile, we assume that their AoDs 
follow the i.i.d. uniform distribution Z/f[0,27r) Q- Due to 
the limited dominant scatterers, the limited AoDs appear the 
SCS. Then, based on ( fTOl i and according to the definition of 
asymptotic orthogonality Q, we have 


proportional to l/V* when x > 0 increases ifTOl . we have 


lim = lim 

M-^oo II J I M^c 

fl 

-Ea||^^pqn|| Jo(am) 


M —1 

|-Ea{^^pgn}| E Jo(am)^ 


< lim- 


SpSq lim 


M-l 

X; JoiamY 

m=0 


(15) 


= 0, 


M^-oo M-s-oo 

where the first inequality is based on Jensen’s inequality. 

Moreover, we investigate the second-order moment of 
^p[TAKq[n]*/M, namely. 


® E ( E E ^ ( 1 -^) 

m=0 Sp=l 

<3 md sin(0o'^) . 

X E 

s„ =1 


( 11 ) 


var { I | _j;| gpWg.W* (jg) 

where var{-} is the variance of a RV. Furthermore, we have 


where p and q denote two arbitrary users with 1 < p < <? < 
K. Based on (HUl, we will derive the first and second-order 
moments of gp[n]gg[n]*/M. Taking the statistical expectation 
of (HUl with respect to AoDs, we obtain 

Ee {gp[n]g9 HVM} 

M — 1 Sp 27r mdsin.(6„^) \-y 

= E E 

m=Q Sp = l 0 

27r md sin(6") 

X E 

Sq=l 0 

Sp Sq M — 1 rj 

= E E yyv&if E Jo(2^^(i-^)) 

Sp = lsg=l m = 0 

M-l 

- E Jo (am)", 


IgpHgqW*! 




n — 0 St) — 1 Sn — 1 


“e" E E 


= 0 3^=1 s' =1 


and we can obtain its expectation for AoDs as follows 


(17) 


M 


-IM-lf Sp 


M 


m=0 


( 12 ) 

where Eg {■} represents the expectation operator for random 
variables (RVs) and subjected to the 

i.i.d. 14 [0, 27r), Jq (•) is the zero-order Bessel function of the 
first kind with the definition Jq (x) = ^ ifTOl . 

Znl=ap’^ PrTp, Znl=aq'’Pu\, « = 27rf (1 - ^), and Vpqn = 
Y.sp=i Ef,Ei cip”f^np{aq”P'nlY■ Besides, due to fc » fe in 
typical massive MIMO systems, we have <C 1 

for 1 < n < and a > 0, where is the carrier frequency. 
According to (fTU . the first-order moment of gp[n]gg[n]*/M 
can be expressed as 


^ | |gpWg^Nl- | g g |4-,|"|4’,|"jo(a(m-m'))^ 

J m=0m'=0l Sp = lSq = l 

+E E E |znp|^( 2 n 9 )* 2 nVo(a(m-m'))Jo(am)Jo(am') 

Sp = lSg = l s' =l,Sg7^Sq 

Sq Sp Sp / / \ * 

+E E E ( 2 np ) 2 npJo(a(m'-m))Jo(am)Jo(am') 

+E E E E 2nyo(am)"jo(am')"l . 

Sp = lSp = l.Sp#SpS5 = ls^ = l,s'/ J 

(18) 

Similar to (fT4l) . we can get the following inequality 


J |gpWgg W| [ _jp]jp\ IgpWggN’ 


M-LM-l, 

<lkE E{jo(a(m-m')) 

m.=Ch-y-i' — fr- 


. . m=0m'=0'' 

+ {Sq+Sp—2)Jo (a(m—m')) Jo (am) Jo {am') 

+ {Sp—l){Sq — 1) Jo(am)"jo(am')"| . 

Similar to (fTSl) and (fTSl) . we can further obtain 


(19) 


fgp[’^]g9['fi]*l _ IP fip fgp[™]g<i['r‘]*ll _ 
\ M j—iSa 11 


Ea{upqp\ E 


(13) 

where E {■} is the expectation operator for all RVs, Ea {;} is 
the expectation operator for RVs and 

and path gains and AoDs are mutually independent. 

Since the large-scale channel fading has no impact on the 
asymptotic orthogonality of channel vectors, without loss of 
generality, we consider channels with normalized energy, i.e., 
Efp=i |Q^?r = l EspEi |ct 9 ‘’|^ = l- Then we have 

|Up9ri| < E E |a?Mnp| 1 ( 09 " Mn,)*! < y/SpSq, (14) 

Sp = 1 Sg = 1 

where the second inequality utilizes the Cauchy-Schwarz 
inequality ifTOll . 

Consequently, based on the discussion above, as well as the 
boundedness (| Jq (x)| < 1 for x > 0) and the delay of Jq (x) 


jta P = jta p. P, =0. 

( 20 ) 

Then we have lim var {gp[n]gq[n]*/M }= 0 by substituting 
(flSl l and (I 20 I 1 into^U. Consequently, the asymptotic orthog¬ 
onality of row vectors of the channel matrix G[n] is proven. 

It should be pointed out that we can further exploit 
the statistical characteristics of path gains to simplify (O 
and (HSJ. Specifically, if we consider the S'fc-path channels 
with following the i.i.d. CAf{0,l/Sk), we have 

= T E^iVpqr,} = 0, £;{zEp} = 0, 
and E { 2 ^} = 0. Therefore, (fTSl l and (fThl l can be further 
expressed as 

M-l 

E {Sp[n]gq[n]*/M}= Y1 Jo{am)^Ea{vpqn}/M ^ 0, (21) 

m=0 

var{gp[n]g5[n]7M}=j:i, |Ee ||gp[n]gJn]*|"/ME| 

M-l M-l - (22) 

= E E Jii{a{m-m')f/M\ 

m = Q m'=0 
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Fig. 1. |E{gp[n]gq[n]*/M} | and var {gp[n]gq[n]*/M} against M and d. 

IV. Numerical Results 

In simulations, the widely adopted ULA at the BS was 
considered, system carrier was /c = 2 GHz, system bandwidth 
was fs = 10 MHz, the OFDM size was N = 4096, 
and the guard interval length was Ng = T^naxfs = 64 to 
combat the maximum channel delay spread of Tmax = 6.4 
m. Besides, we adopted Rayleigh fading sparse multipath 
channels with typical S = 6 paths model a, where the 
normalized maximum delay spread is Ng = 64 and path gains 
follows the i.i.d. CAf{0,l/S). While for the ideal Gaussian 
channel matrix, the large-scale channel fading is removed 
(D = Ifc), and elements of small-scale channel fading matrix 
H follow the i.i.d. CAf{0, 1). 

Fig. m shows the simulated and analytical 
\E {gp[n\gq[n\* /M} I and var {gp[n]g,[n] against M 

and d. Clearly, the analytical (ImTi and (l2^ have the excellent 
tightness with the simulation results, which confirms the 
validity of our theoretical analysis on the first and second- 
order moments of the inner product of channel vectors. 
Moreover, var{gp[n]gq[n]*/M} decreases when d increases, 
which implies that the larger d can lead to the smaller 
inter-user interference. Besides, both analytical and simulated 
curves confirm the asymptotic orthogonality of different 
channel vectors over sparse multipath channels. 

Fig. 12] provides the cumulative distribution function (CDF) 
of the smallest/largest eigenvalue of GG* for small-scale 
MIMO systems (M = K = 6) and massive MIMO systems 
(M = 128, K = 6). Meanwhile, we plot the CDF of the 
smallest/largest eigenvalue of GG* under the ideal Gaussian 
channel matrix assumption as the benchmark for comparison. 
Clearly, the gap of eigenvalue distributions of GG* over 
sparse multipath channels and that under the ideal Gaussian 
channel matrix is negligible, and this gap is narrowed when 
d increases. Besides, the simulated CDFs of eigenvalue dis¬ 
tributions of GG* over sparse multipath channels conform to 
the experimental results provided by Fig. 6 of m. 

Fig. [3] shows the capacity of massive MIMO systems with 
M=128 and A=16 over sparse multipath channels and that under 
the ideal Gaussian channel matrix assumption. The capacity 
gap between massive MIMO systems over realistic sparse 
multipath channels and that under the ideal Gaussian channel 
matrix assumption is small, and this gap is narrowed when 
d becomes large, since the inter-user interference is reduced 
when d becomes large as confirmed in Fig.[T] Note that we do 
not consider the antenna coupling, which may further degrade 
the attainable capacity when d is small. The small capacity 


Eig. 2. CDEs of smallest/largest eigenvalue of GG* in MIMO systems. 
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Eig. 3. Comparison of capacity per user (bits/channel use) against the antenna 
spacing A in the massive MIMO system with M = 128 and K = 16. 

loss conforms to the experimental results that the capacity of 
implemented massive MIMO always approaches the capacity 
under favorable propagation condition IT]. Besides, the small 
capacity loss implies that we may realize BS antenna array 
with small form factor (e.g., d < 0.5A) at the cost of a small 
capacity loss if the antenna coupling can be settled well. 

V. Conclusions 

This paper has investigated the propagation condition and 
system capacity of massive MIMO over realistic sparse chan¬ 
nels with SCS in both time domain and angle domain. We have 
proven the equivalence between the favorable propagation 
condition desired by massive MIMO and realistic channels 
with SCS. Specifically, we have derived the first and second- 
order moments of the inner product of different channel vec¬ 
tors over sparse multipath channels. Moreover, the asymptotic 
orthogonality of channel vectors over massive channels with 
SCS is proven. This proof indicates that the capacity of 
massive MIMO over realistic massive MIMO channels with 
SCS can approach that of the ideal Gaussian channel matrix 
assumption when the number of BS antennas is unlimited. 
The asymptotic orthogonality of channel vectors has been 
verified by comparing the simulated and analytical first and 
second-order moments of the inner product of channel vectors. 
Moreover, the simulated condition number of channel matrix 
and achievable capacity of massive MIMO systems over sparse 
MIMO channels are consistent with experimental results. The 
negligible capacity gap between massive MIMO over realistic 
sparse channels and that under the ideal Gaussian channel 
matrix assumption may inspire us to realize massive MIMO 
with better performance by exploiting the channel sparsity. 
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